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Abstract. We present theoretical schemes for performing high-fidehty one- and two- 
qubit pulsed gates for a superconducting flux qubit. The "IBM qubit" consists of three 
Josephson junctions, three loops, and a superconducting transmission line. Assuming a 
fixed inductive qubit-qubit coupling, we show that the effective qubit-qubit interaction 
is tunable by changing the applied fiuxes, and can be made negligible, allowing one to 
perform high fidelity single qubit gates. Our schemes are tailored to alleviate errors 
due to 1/f noise; we find gates with only 1% loss of fidelity due to this source, for 
pulse times in the range of 20-30ns for one-qubit gates (Z rotations, Hadamard), and 
60ns for a two-qubit gate (controlled-Z). Our relaxation and dephasing time estimates 
indicate a comparable loss of fidelity from this source. The control of leakage plays 
an important role in the design of our shaped pulses, preventing shorter pulse times. 
However, we have found that imprecision in the control of the quantum phase plays 
the major role in the limitation of the fidelity of our gates. 
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1. Introduction 

Superconducting circuits containing Josephson junctions [1^ 2j are widely recognized to 
be promising systems for the physical implementation of quantum bits. These systems 
can be made and operated using well-established experimental techniques, and they 
have the clear potential, in principle, for scalability. Important experimental milestones 
in coupled superconducting qubit s, including the observation of two-qubit gates [3^ 4j 
and the measurement of entanglement [5j, have already been reached. But because 
superconducting qubits are condensed matter systems, it has been a hard task to isolate 
them from their environment. Because of that, these systems tend to suffer from short 
coherence times, which has imposed serious limitations on achieving very high fidelity 
gates. 

In this paper, we report a theoretical study of a universal set of one- and two-qubit 
gates implemented using only shaped dc flux pulses for an oscillator stabilized flux qubit. 
We introduce a simplified but accurate model to describe the dynamics of the lowest 
states of the qubit as a function of the external control parameters. This model provides 
a simpler way to analyze the physics of the problem, helping to make it easy to see what 
operations are necessary for performing the desired quantum gates. As we shall see, 
there is a trade-off between the speed of a gate and the amount of leakage produced by 
it (leakage = evolution of states out of the 0-1 computational basis). Smart choices for 
the shape of the pulses are required to keep leakage at a tolerable level, such that the 
gate fidelities are not compromised by this process P [7] . We have been able to keep 
leakage at the 0.12% level for gate times of order of few tens of ns. 

The other important goal of the present work is to consider the effect of low- 
frequency noise during the gate operation and in the memory state. When designing 
quantum gates, we search for the best pulse paths such that the loss of fidelity due to 
the low- frequency noise is minimized. We find that for all gates of interest, the loss of 
fidelity due to low frequency noise is never greater than 1%. Achieving this low level of 
infidelity requires a careful use of symmetries (so that errors accumulated in the first half 
of a pulse can be cancelled out in the second, for example) and of "sweet spots" ^ Q\ 
(points in the control parameter space that are first-order insensitive to fluctuations). 

The last part of this paper analyzes a pulsed two-qubit gate, the controlled phase 
gate. An interesting feature of the coupling used in this gate scheme is that, although 
the physical inductive coupling between the qubits is assumed fixed, the effective qubit- 
qubit interaction is tunable as a function of the control parameters on both qubits, due 
to the change of character of the bare qubit states as the control parameters are varied. 
In fact, the effective qubit-qubit interaction can be made negligible in large regions of 
the flux space, allowing us to perform high fidelity one-qubit gates for the two-qubit 
system without extremely stringent control over electrical cross-talk. 

The outline of this paper is as follows. In sections [2] and [3] we discuss the physics of 
our system, and present the simplified model used to simulate the dynamics of the lowest 
lying levels of an oscillator stabilized flux qubit. The Appendix carries out a detailed 
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derivation of the system Hamiltonian, and the regime of vahdity of the simphfied model 
is discussed. Section [4] presents the different schemes designed to perform the foUowing 
basic single-qubit operations: measurement in the standard basis (0/1), measurement in 
the conjugate basis (+/-), Z-rotation gates, and the Hadamard gate. The fidehty of aU 
these gates is given as a function of unwanted shifts from the optimal point of operation 
in flux and time synchronization of the pulses. In addition, we characterize the noise 
through the operator-sum representation of the system superoperator. In section [5| 
the two-qubit system is analyzed. The form of and the reasons for an eflFective tunable 
interaction are discussed. Then, a gate in the equivalence class of the controUed-Z gate is 
proposed, and its fldelity and a characterization of the nature of the noise are presented. 
Finally, section [6] gives some conclusions. 

2. System Hamiltonian 

Our analysis will be focused on the IBM qubit [TOl |TT]. This device, shown in Figjl] 
consists of a bare qubit^ which is a type of flux qubit containing three Josephson junctions 
and three loops, and a high-quality superconducting transmission line [121 [HI E]- 
The bare qubit is subject to external control via flux lines which change the total 
magnetic fluxes threading the loops. As previously reported p^, and summarized in 
the Appendix, the bare qubit has a gradiometric structure, and hence its behavior, 
to good approximation (see Appendix for the validity of this assumption), is only a 
function of the diflFerence of the magnetic flux in the two large loops, which we denote 
as e = $ — $p. Whenever e is an integer multiple of the flux quantum $o = h/2e^ the 
system potential has a perfectly symmetric structure. These lines in the flux space are 
referred as the "S Lines" (S for symmetric). 

The bare qubit has two control parameters that deflne its flux space: one, e, causes 
departures from the parity symmetry just described, and the other, the so called control 
flux changes the structure of the qubit potential, varying the height of the potential 
barrier between the two classically stable states. As $c is changed, the system potential 
passes from a double- well to a single- well structure. This change in the form of the 
potential as a function of a control parameter produces an essential feature in the 
quantum behavior of the system: a well-deflned change in the character of the qubit 
eigenstates. While for the double-well regime the qubit eigenstates are almost localized 
orbitals in the left and right wells, as the single well regime is approached the qubit 
eigenbasis goes over to delocalized states, being close to states that are symmetric and 
antisymmetric with respect to reflection around the midpoint of the potential. 

This change of the character of the states happens over a very small interval of 
because of the exponential increase of the amplitude of tunneling between the wells as 
the barrier between the two wells is decreased. This transition interval is referred as the 
''portal" [ini [H], and pulsing the flux parameters through this portal can create some 
of the elementary actions needed for the construction of quantum gates. For example, 
a non-adiabatic pulse through this region can create superpositions between states |0) 
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Figure 1. (a) Schematic layout of the hare IBM qubit. This device consists of three 
loops and three Josephson junctions (indicated by red crosses), and operates in the 
flux regime, (b) Picture of the qubit. The white tracks highlight the bare qubit (a), 
and the transmission line coupled to it (white tracks leaving the right-hand side of 
the picture). In order to operate the qubit, flux lines (highlighted in black tracks) 
are used to change the total fluxes threading the loops. Readout SQUIDs (two such 
structures are shown at the left-hand side of the picture) perform the measurement 
of the state of the qubit. (c) Proposed scheme of the two qubit system. The qubit- 
qubit interaction is assumed to arise via the indicated mutual inductance between the 
two big loops [14 . This results in a qubit-qubit interaction Hamiltonian of the form 
CFz ® CFz. Physical parameters for this qubit: the capacitances and critical currents of 
the junctions are assumed C =10fF and Ic = 1.3/iA; Lt = 5.6nH, Li = 32pH and 
Ls = 680pII are the transmission line, the small and big loop inductances, respectively. 
The mutual inductance between qubit-transmission line, small loop-control flux line 
and big loop-bias flux line are respectively 200pH, O.SpH and 0.5pH. The transmission 
line is designed to have a fundamental mode frequency of cot = 27r x 3.1GHz. Finally, 
the qubit-qubit mutual inductance is M =12pH. 

and |1). 

The fundamental mode of the open-ended transmission line acts as a harmonic 
oscillator of frequency ujt coupled to the bare qubit. The presence of that structure 
modifies the quantum behavior of the qubit when the energy splitting of the ground 
and first excited state of the bare qubit is comparable to or larger than Hut. In this 
regime, the two lowest-lying states of the system both have the bare qubit in its ground 
state; they differ only in their transmission line quantum number. By tuning the energy 
splitting of the ground and first excited state, one can move information stored in the 
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bare qubit to the transmission line and vice- versa. When the qubit has been transferred 
to this transmission-hne embodiment, we say that it is parked. Parking results in a very 
useful stabilization of the 0-1 frequency as a function of changes in {$C5 In addition, 
the quantum coherence times while parking are seen to reach several /isfUj. Thus, the 
parking regime will be used as the memory state of the qubit; it will stay in this state 
when it is awaiting operation. Far away from parking, when the 0-1 energy gap is much 
smaller than Iiwt^ the transmission line does not play any role in the dynamics of the 
0-1 states, which are just those of the bare qubit. Since the coherence times here are 
expected to be of the order of only tens of ns, this regime should be avoided; this part 
of parameter space will only be used for state measurement. 

So far, we have given a qualitative description of the qubit dynamics. This 
description is made quantitative with the methodology introduced by Burkard, Koch 
and DiVincenzo (BKD)p^. Using network graph theory, BKD developed a universal 
method for analyzing any lumped element electrical circuit containing Josephson 
junctions. The result of this theory is a mapping of the circuit dynamics to that of a 
massive particle in a potential, whose mass tensor and degrees of freedom are associated 
with the system capacitances. The system Hamiltonian in this formulation is given by 

ns{t) =lQ5c-iQe+(^^) t/(<^,t), (1) 

1 2lT 

u{ip, t)= -J2 ^ + 2^"^^°^ + w^"^^^ * ^^^^^^ + * 

i 

Here Lj^ = (^^^ h-ii where Ic-i is the critical current of junction i. The diagonal matrix 
C contains the capacitances of the system. The topology of the circuit is encoded in the 
matrices Mq, N and S (see Appendix for a more detailed presentation of BKD theory). 
The first term of the potential Eq. ^ represents the energy due to the presence of the 
Josephson junctions. The second term is associated with the inductive energies of each 
branch of the circuit. The last two terms take into account coupling to external sources 
of magnetic fiux and current sources I^. Quantization of the system is introduced 
by imposing the canonical commutation relation for the variables of charge, Qc, and 
phase, (p: 



ihSij. (3) 



The analysis of the system potential Eq. ^ for our qubit reveals that, instead of 
working with the real applied fiuxes = {^c, ^p}^ it is more convenient to introduce 
''non-orthogonal" fiux coordinates e = $ - $p and $c = + 2{l^I'm^-Mi^) + ^v) 
(where Li stands for the small loop inductance, L3 to big loop inductance, and M15 and 
M35 to the small-big loop and two big loops mutual inductances, respectively). In fact, 
as presented in the Appendix, the potential has a definite symmetry as a function of 
{^c)^}) rather than {$C5^}- Since the potential symmetry is a key feature of system, 
and we will explore it when designing our gates, an accurate calculation must consider 
that fact, even though the term proportional to 2{L^l'M^-\iir,) ^ 0.024 represents a small 
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correction to the real applied flux $c- From here onwards, we will only refer to the pair 
of effective fluxes {$C5 

Since our qubit has four capacitances (three associated with the Josephson 
junctions, and one representing the fundamental mode of the transmission line), BKD 
theory leads us to a four-dimensional potential, for which direct calculations and 
analysis are more diflicult. Thus, we follow the procedure of [13] to reduce the system 
dimensionality to two, one representing the bare qubit and the other representing the 
transmission line. The procedure involves organizing the degrees of freedom into "fast" 
coordinates, in which the potential rises very steeply (such that the system dynamics 
are frozen into ground state along these direction), and ones that are "slow". Then, 
using a Born-Oppenheimer approach, the fast coordinates are traced out, resulting in 
small modiflcations to the remaining slow-coordinate potential energy. 

Fig. [2] shows the flrst seven levels of the system spectrum, calculated following the 
above described steps, as a function of for three different values of the bias flux 
6. Note that, by convention, the ground level is always at = 0. It is clear that the 
equally spaced harmonic oscillator levels cut through this spectrum for all values of $c 
and e; since the transmission line sees the control flux only via interaction with the bare 
qubit, its energies are very stable except in the vicinity of energy crossings. Observe 
that the lowest-lying states in the "parking" regime, at high values of involve only 
transmission line quantum numbers (i.e., they are the states of a harmonic oscillator). 
This is true because the energy splitting of the eigenstates of the bare qubit becomes 
much larger than the hujT energy splitting of the transmission line states. 

In addition, because of the interaction between the bare qubit and the transmission 
line, another structure present for each e is an observed avoided crossing gap between 
the bare qubit and the transmission line states, occurring close to $c = 1.45$o- As one 
can see, the bias flux e plays an important role for small values of where the system 
potential has a double- well structure, and the bare qubit energy splitting is smaller than 
IijjJT' The flrst plot presents the case on the S line, 6 = 0. Since at the S line the system 
potential is symmetric, the qubit states are symmetric and antisymmetric superpositions 
of the degenerate localized orbitals of the left and right wells of the potential. As we 
move away from the S line (e 7^ 0), the symmetry of the potential is broken and the 
localized left and right orbitals are no longer degenerate. This explains the appearance 
of the gap between the ground and flrst excited states in the second and third plots 
at small values of $c- One important feature of the system potential is a symmetry 
in going from +e to —e (see Appendix); under this transformation, the energies of the 
left and right states are interchanged. Finally, in the third plot it can be seen that the 
flrst excited state stabilizes at the transmission line frequency for $c < 1.43$o- This 
occurs because for most values of $c the bias term in the Hamiltonian is large enough 
by itself to make the bare qubit energy splitting larger than fwoT (but not much larger) . 
As a result, for small $c the flrst excited state corresponds to one excitation of the 
transmission line, and the second excited state to the excited state of the bare qubit. 
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Figure 2. The first seven system eigenlevels, as a function of the control flux <l>c, 
for three different values of the applied flux e. The qubit is encoded using the the 
lowest two eigenstates. For large values of the energy splitting of the bare qubit 
states is much larger than Hujt^ so that the lowest lying qubit states have purely 
harmonic oscillator character, (a) The level structure on the S line. In this case, the 
system potential has a symmetric double- well structure for small values of control flux, 
consequently the ground and first excited states are nearly degenerate there, (b) and 
(c) present cases away from the S line. Here, the states |L) and \R) are no longer 
degenerate and the system has a gap between the ground and first excited states for 
small values of control flux. 



3. 4- level model 

Although it would be possible, with considerable computational effort, to design 
and simulate the desired quantum gates by direct evaluation of the time dependent 
Schrodinger equation for the circuit Hamiltonian Eq. Q, a considerable economy 
is achieved by introducing a simplified model that correctly describes the lowest 
eigenstates, since we will be only interested in their dynamics for our quantum gates. 
We start the derivation of such a simplified model by using first order perturbation 
theory to treat the bare qubit coupling with the flux line and with its transmission line. 
The resulting Hamiltonian due to that approximation has the form of a biased two-level 
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Figure 3. Coefficients of the one-qubit Hamiltonian Eq. Q as a function of the 
control flux The solid curve gives the tunneling amplitude between the ground 
and first excited states of the system. The dashed curve gives the coefficient of the 
coupling between the bare qubit and the transmission line, and dot-dashed curve gives 
the bias term. For small values of ^c-, the tunneling amplitude becomes exponentially 
suppressed due to the very large barrier between the two wells of the potential. As one 
increases <l>c, the potential barrier decreases and the positions of the minima become 
closer together. Around = 1.45<l>o the barrier vanishes entirely; from this point 
onward, the potential has a single well, as sketched in Fig. [2] Inset: detail of the 
tunneling amplitude for the single- well regime. 



system coupled to a harmonic oscillator (see Appendix for a detailed derivation) 

H = -^A($,)a, + + hwra^a + g{^c){a + d^)a,, (4) 

where are the Pauli matrices, a^lL) = — ^z\R) = +1^)^ {a'^ ^ a} are canonical 
bosonic creation and annihilation operators. The states \L) and \R) represent the 
localized orbit als found in the left and right potential wells. 

Because of the change of the state character as a function of the amplitude 
of tunneling A, the bias term coefficient b and the qubit-transmission line coupling g 
also become control flux dependent. Fig. [3] presents their behavior as a function of 
$c- As one can see, the tunneling amplitude A becomes negligible for small values of 
$c- This happens because, in this regime, the very large barrier between the two wells 

lOOGHz) in the double- well potential exponentially suppresses the tunneling between 
their lowest states. However, as the value of $c increases, the barrier decreases and the 
two minima become closer. Consequently, the amplitude of tunneling rapidly increases, 
and the left and right orbitals start to become more and more delocalized. Around 
$c ^ 1.45$o the barrier vanishes rather abruptly and the bare qubit potential enters 
a single- well regime. The region around $c ^ 1.45$o deflnes the portal of the system. 
The dashed and dot-dashed curves present the coefficient of the bare qubit-transmission 
line coupling g^ Eq. (A. 20), and the bias term 6, Eq. (A. 19), respectively. Once the bare 
qubit reaches the single well regime, those coefficients become less sensitive to changes 
in since the nature of the states barely changes as one passes through this regime. In 
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Figure 4. The first system eigenstates as a function of control flux obtained 
through the Schrodinger equation solution of the exact potential Eq. ([2| (solid curve), a 
simple truncation of Hamiltonian Eq. Q at the first excited state of harmonic oscillator 
(circle-symbolized curve), and the solution of the 4-level propos ed, Eq. ([5| , using the 
shifts 8 = —g/hwT (triangle-symbolized curve) and s = —g/\/u;^-\- (dashed red 
curve), Eq. (|6|, for three different values of the bias flux e. The simple truncation fails 
when g becomes appreciable (small values of ^c)- In this regime, the shifted harmonic 
oscillator states are the preferred system representation. The 4-level model using the 
shift 8 = —g/hujT does not provide the parking harmonic stabilization, since it does 
not consider the change of character of the bare qubit. At last, the 4-level model using 
the conditional shift Eq. (|6| gives a very fair approximation of the lowest-lying levels 
for all regimes of the bias and control flux. 

addition, the inset shows the amplitude of tunneling can rapidly reach values of several 
tens of GHz in the parking regime, which are, in general, much larger than the typical 
value of the transmission line energy splitting Hujt^ hence the lowest-lying states in the 
parking regime involve only transmission line states. 

Unfortunately, even though the Hamiltonian Eq. ^ already represents an 
important simplification for the system description, since it has reduced the bare qubit 
Hilbert space to that of a two-state system, we cannot obtain an analytical solution 
for its eigenstates/eigenvalues. Thus, we still have to deal with an infinite set of states 
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due to the harmonic osciUator. As we shaU see, for the purposes of our work, we must 
have not only the correct dynamics of the states |0) and but also an exceUent 
agreement for the minimal gap between the computational basis {|0), |1)} and the rest 
of spectrum of the system. Because of that, the truncation of the harmonic oscillator 
Hilbert space at its first excited state does not work as a fair approximation of our 
system Hamiltonian. As presented in Fig. [4| the simple truncation of the harmonic 
oscillator Hilbert space fails to give a good description of the gap between the |1) and 
1 2) states. This happens because at the regime of small values of $c the function g 
is appreciable. Hence, the shift applied to the oscillator due to the interaction with 
the bare qubit becomes important. Consequently, the most adequate description of the 
system is obtained using the representation of shifted harmonic oscillator states, where 
the truncation of the harmonic oscillator Hilbert space should give better results. 

In fact, a more careful inspection of the Hamiltonian Eq. ^ reveals that its 
harmonic oscillator part has a form of a well-known shifted oscillator. Thereby, changing 
to the representation of the shifted harmonic oscillator states should give us a better 
picture of the system dynamics, and the appropriate basis in order to perform further 
truncations of the Hilbert space. Nevertheless, the shift applied to the harmonic 
oscillator is dependent on the bare qubit state, as can be seen by the term g{^c){ci+d))az 
of Eq. This important feature of the system leads us to introduce the unitary 

transformation D the conditional displacement operator of the 

system. Observe that the operator D does not commute with the spin part of the 
Hamiltonian Eq. Q. Thus, changing the Hamiltonian representation to the shifted 
harmonic oscillator states, H HD^ and then performing the truncation of the 

harmonic oscillator Hilbert space at its first excited state, we arrive with the following 
simplified 4-level model 

+ Hujtc^c + (T^c(^($c) + s^Hujt) + c^a^(^($c) + sHujt), (5) 

where we have defined the operators c = |0ho)(1ho| and = |1ho)(0hoL which |0ho) 
and |1ho) represent the ground and first excited harmonic oscillator states, respectively. 
During the procedure described above, we have introduced an ad hoc parameter, the 
shift parameter 5, which can be used to parametrize the Hamiltonian Eq. ^ in order 
to obtain the closest level dynamics as possible to that determined by the Hamiltonian 
Eq. Q. 

A natural choice for the parameter s would be the standard value s = —g{^c)/^T-> 
which leads to the cancellation of the last two terms of Eq. ([5]), and the Hamiltonian 
diagonalization when A 0. However, since the shift imposed on the harmonic oscillator 
is conditioned on that of the bare qubit state, and the form chosen for s does not take 
into account the change of character of the bare qubit state, one should expect it would 
fail when the regime of high tunneling amplitude is reached. Indeed, as can be seen 
in Fig. [4| this choice for the parameter s only gives the correct description when the 
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tunneling amplitude is negligible. In fact, the analytical solution of Eq. ^ using the 
above shift reveals that, at the limit A ^ oc, one should expect cjqi oc, what is in 
complete opposition with the parking stabilization expected for the system eigenstates. 

Therefore, a smart choice of the parameter s has to consider the fact of the harmonic 
oscillator shift is conditioned to the bare qubit state, and that its state character changes 
as a function of $c- As we already know, at the limit A ^ 0, the parameter s should 
asymptotically reach the value s —> —g{^c)/^T) since it correctly decouples the spin 
and harmonic oscillator degrees of freedom in the Hamiltonian Eq. ([5]). In the other 
limit, A ^ oc, one should expect observing no shifts to be imposed on the harmonic 
oscillator, since the system states would be frozen at the bare qubit ground state. Thus, 
we expect a reasonable interpolation for the parameter s between the two regimes is 

si^c) = - r^^i^ ,' (6) 

In fact, as presented by the dashed curve in Fig. [ij the form Eq. ^ for the 
parameter s gives a very fair approximation for the lowest level dynamics, in particular 
the energy splittings cjqi and (jJ^. In addition, from the analytical solution of Eq. ^ 
using Eq. ([6]), we can check that the cjqi has the correct harmonic oscillator stabilization, 
z. e. A ^ oc ^ Uqi ujt- 

It is worth pointing out that, although the model Eq. ^ does not give the correct 
parking stabilization for the system second excited state |2), since this state involves 
the second excited state of the transmission line, for the purposes of our work, it turns 
out that is not a limitation for the model. Indeed, the lack of parking stabilization for 
the second excited state |2) would lead to underestimations of the transitions between 
the state |1) and |2), and consequently wrong leakage estimations. However, as we 
are going to assume ujt of the order of several GHzs, the energy splitting uji2 for the 
parking regime is big enough to suppress those transitions during our dc shaped pulse 
operations. Furthermore, since at the parking regime the ground state, and the first and 
second excited states have the same nature [i. e. they are the states of the harmonic 
oscillator), observing the induced 0-1 transitions due to the dc pulse operations also 
gives a measurement of the expected 1-2 transitions in this regime. Thus, we end with a 
very controllable four dimensional Hilbert space model, spanned by the ground and first 
excited states of the bare qubit and transmission line, that accurately mimics the main 
features of the system dynamics determined by the exact circuit Hamiltonian model Eq. 
Q. From now onwards, we will only use the Hamiltonian Eq. ([5]), with Eq. ([6]), when 
designing and simulating the quantum gates. 

Completing the system description. Fig. [5] presents the plots for the frequency 
difference between the ground and first excited state, cjqi, and the frequency difference 
between the first two excited states, 6^12, as a function of {^c, e}- As mentioned before, 
the plot quantifies the minimal gap between the computational basis {|0), |1)} and 
the rest of spectrum of the system, giving a measure of how likely the system is to 
undergo leakage during the sweep of a pulsed gate. The cjqi plot shows how fast the 
relative phase of the two computational basis states advances if the system is held at 
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some flux values. We envision that there wiU be a very stable master clock at the 
transmission-line frequency, and that errors in the accumulated phase diflFerence are 
unlikely if the system is held in parking, where its phase advance is synchronous with 
this clock. But when cc^oi departs far from cj^, phase accumulation with respect to 
the reference is very fast, and we expect (and our numerical studies confirm) that 
the system is much more susceptible to phase errors in this regime. These plots 
provide a good "map" for designing the one-qubit gates, since they indicate the regimes 
where one should expect an appreciable amount of leakage and very fast relative phase 
accumulation. 

4. One-qubit Gates 

Before we move on to the discussion of each gate individually, it is worth summarizing 
some features present in all of them (including the two-qubit gate). First the memory 
and measurement points must be defined. The parking point is taken to be at 
$c = l-6$0 5 the measurement point to be at $c = 1.4$o- These choices are 
optimal, in the following sense: For the parking point, the important feature to be 
considered is how much the cjqi frequency chang es as a function of external controls. 
Since we envision the qubit spending long times at this position, it is imperative that 
the frequency deviation ^c^Jqi have very small values for reasonable fiux shifts due to 
the noise. We calculate that the frequency sensitivity is 8ujqi/5^c ^ 10MHz/$o at 
the chosen memory point, which will give acceptably small memory error. In addition, 
inaccuracies in the transmission line fundamental model frequency can be corrected 
using a vr-pulse stabilization scheme [17j. This operation (not analyzed in detail here) 
consists in initially applying a tt rotation gate; the system is then left to evolve in a 
free evolution for the same amount of time previously spent at the memory state, and 
finally a new tt rotation is applied. As a result, undesired Z-phase accumulated at the 
parking point ends as a global phase of the system state. 

For the measurement point, we must make sure the left and right states are 
experimentally distinguishable, and that, since the potential has a periodic structure 
[13j, the potential barrier between different pairs of minima is high enough to avoid 
tunneling between their states. At $c = 1.4$o we found the barrier between the two 
principal minima is of the order of lOTHz, and the barrier to other minima is even 
higher (c. 20THz). 

Another common feature for all of our shaped pulses is our design methodology. 
We shape our dc pulses using simple sums of tanh functions. As we shall see, our gates 
can be divided in several parts, and with each part is associated a function of the form 
5 tanh ( ^~^^^^ ), where 5^ r and tr^ax determine the fiux excursion, the maximum rate 
and its time position, respectively. Thus, with those parameters, one can adjust the 
rate and the fiux excursion of each part of the gate, in order to optimize the gate. It 
turns out that our designed gates do not require maximum fiux slew rate and bandwidth 
higher than 7 x 10^$o/s and IGHz, respectively. 
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Figure 5. Plots of the energy splittings (a) UO12 and (b) cjqi for the flux space 
{<l>c,e}. The plots provide a good "map" for designing the one-qubit gates, since 
they indicate the regimes where one should expect an appreciable amount of leakage 
and very fast phase accumulation, (a) The UO12 plot quantifies the minimal gap between 
the computational basis {|0), |1)} and the rest of spectrum of the system, giving a limit 
on the rate at which the system can evolve without appreciable leakage, (b) The cjqi 
plot indicates the phase accumulation rate with respect to the reference. In regions 
of very fast rates, a very precise control of the external applied fiuxes is required in 
order to avoid phase noise. Illustrated with dots is the path in the fiux space used 
to implement the Hadamard gate. The number of dots indicates the time spent when 
passing through that region. 



We measure the gates' fidelity using the entanglement fidelity [ 18] 

where = Ul^^^iUreai{5^c,5e,5t), with Uideai and Ureai{5^c,5e,5t) representing the 
ideal gate and the final achieved transformation, respectively, is an equal distribution 
of the computational basis states |0) and |1): = ||0)(0| + ||1)(1|. Finally, we estimate 
the probability of leakage using the projection of an arbitrary evolved state, |^^/) = 
U\^i), in the 0-1 computational basis: 1 - (^/|(Po + = 1 - |(i^o + Pi)U\^i)\^, 

where Pj represents the projection operator in the state 
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Figures |6| [7] and [8] show the proposed measurement gates, phase gate, and the 
Hadamard gate, whose matrix representation in the {|0),|1)} basis are respectively 
given by 

/ 1 n \ 

(8) 
(9) 



U^phase - I Q ^io. 1 ' (10) 
UHadamard ^ _^ I • fill 

Where ^oi and are arbitrary phases, whose values are not relevant for the gate 
implementation, provided that the gates Eqs. ([s]) and ^ are followed by a projection 
into the basis {|0), |1)}. The relative phase 9z is the parameter that defines the phase 
gate, Eq. (10). Finally, since the states |0) and e^^°|0), and |1) and e^^^ll) are physically 



identical, the convention for the global phases Oq and 9i is chosen such that the physical 



implementation of the Hadamard gate has the matrix representation given by Eq. (11). 
Once fixed, the phase convention is maintained for all other gate implementations. 

In addition, figures |6| [7| and |8] also present the gate fidelities as a function of 
unwanted constant shifts from the optimal point of operation in the whole pulse profile. 
We model the eflFect of 1/f fiux noise by an ensemble of random constant shifts of this 
sort. Our assumption of constancy is accurate for those components of the 1/f noise 
at frequencies below the inverse of the gate time, around lOOMHz. 1/f noise at this 
frequency and higher is not accurately modelled in this way, but at these frequencies 
we believe that other sources of high frequency noise (that is, white noise from the 
resistances in our circuit) become more important than 1/f noise. In our work these 
noise sources are modelled separately using a quantum bath; see Ref. [13j for these 
calculations. 

Our estimates for Ti and T2 [16j times indicate that we can expect a very long 
coherence time in parking, 0(1^), and a very short dephasing time, 10ns, at the 
measurement position. In the portal region, we calculate coherence times of the order 
of hundreds of fis. As we shall see, except in the measurement processes, our gates are 
designed to not go lower than the upper limit of the portal. Thus, in principle, for gates 
of duration times of 20-30ns, the fidelity should not be compromised by decoherence 
much more strongly than by the imperfections explored in these figures: unwanted 
shifts in the applied fiuxes due to low- frequency noise, and temporal shifts between $c 
and e pulses. 

Those shifts are assumed uncorrelated and constant during a single gate operation, 
but random from one "shot" to another. We model the noise as a normalized gaussian 
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Table 1. The expected gate fidelities, as a function of unwanted shifts in 5^ci 
and 5t^ calculated considering each channel noise separately. In addition, the expected 
total net fidelity evaluated in the presence all channels is presented. 



One-qubit system 


Gate 






5e 


5t 


All 


0/1 Measurement 


(Fidelity) (%) 


99.99 


99.99 




99.99 


+/- Measurement 


(Fidelity) (%) 


99.81 


99.84 




99.8 


Phase Gate 


(Fidelity) (%) 


99.999 


99.999 




99.999 


Hadamard 


(Fidelity) (%) 


99.87 


99.47 


99.63 


99.46 



probability distribution, jji{x) = -^e~^^/^^^, which we assume having (at IHz) 6/x$o 
(flux shifts) and 6ps (time shifts) as its root-mean-square deviation, a, such that 
approximately 90% of the distribution is found between the values ±10/x$o (flux shifts) 
and zblOps (time shifts); we believe that this quality of control will be readily achievable 
in the lab in the coming years. Indeed, recent experimental results indicate that 
our assumption of having 6/x$o as root-mean-square deviation of the 1/f flux noise at 
IHz is already achievable for Josephson qubit s. 

We note that by using the 1/f rms amplitude at IHz in this model, we are implicitly 
assuming that the 1/f noise is cut off below this frequency. In fact, the 1/f spectrum 
goes much lower (at least three orders of magnitude lower in Ref. [19], and in other 
similar previous studies); but we may assume that in quantum computer operation, 
qubits are frequently taken off-line and recalibrated. Assuming that this recalibration 
takes place once per second gives the IHz cutoff. Our assumption of Gaussian statistics 
merely embodies the expectation that the 1 /f noise arises as the summed effect of many 
independent fluctuators; this is also borne out by the traces of Ref. [19j. 

Table [l] summarizes the average fldelity obtained for the gates discussed considering 
each channel noise separately: (F) = J dSxfi{Sx)F{Sx)^ where F is given by Eq. 
As one can see, the minimal expected fldelity found was 99.47% for the Hadamard gate. 

For the phase, Hadamard and controUed-Z gates, we present the operator-sum 
representation |20j of the system superoperator 

S{p) = J d5xiJ.{5x)Ureai{Sx)pUl^^i{5x) 

^ Y^E.pEl (12) 

k 

where the operators {Ej.} are the operation elements for the quantum operation 8. 
Having the set of operators {E^} permits the characterization of the noise present during 
the system evolution. We flnd that the model noise considered in this work indicates 
that the effective noise for the physical implementation of our universal set of quantum 
gates is heavily biased, i. e., we have found that the effect of phase noise is at least one 
order of magnitude higher than bit-flip noise and leakage processes. New strategies for 
fault tolerant computation have recently been worked out which take advantage of this 
sort of biassing in the noise |2T] . 
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4.1. Measurement gates 

We present two distinct measurement gates, both very useful for effective quantum 
error correction and universal quantum computation. As one can see at their matrix 
representations Eqs. ([s]) and ([9]), by measurement gates we mean the unitary operations 
that are used to prepare the state for the final projection process (measurement). The 
first gate is a measurement in the standard 0/1 basis. That is, starting from parking, we 
are to distinguish whether the system is in its ground state (0 quanta in the transmission 
line mode) or its first excited state (1 quantum in the transmission line mode). This 
gate, shown as the first plot of Fig. [6| works by performing an adiabatic evolution of 
the states |0) and |1) from the parking point to the measurement point, maintaining 
the bias fiux e at a constant value off the S line (e 7^ 0). So, through this adiabatic 
transformation, the qubit states evolve from the configuration IS', 0) (ground state) and 
\S^1) (first excited state), to |I/,0) and 0) (where the first label corresponds to 
the bare qubit states, and the second to those of the transmission line). The matrix 
representation Eq. ([s]) shows that we end with an irrelevant relative phase ^01 between 
the |0) and |1) states. The sources of loss of fidelity are leakage at the avoided crossing 
gap, and 0/1 transitions at the portal. For the path shown, the probabilities of observing 
leakage and 0/1 transitions are 3 x 10~^% and 6 x 10~^%, respectively. As a result, the 
expected total net gate fidelity was found: J dS^cd^e ii{5^c)l^{^^)F{5^c)5^) ^ 99.99%. 

The second measurement gate is in the conjugate basis. It evolves equal 
superpositions of |0) and |1) states at the parking point, |±) = ;^(|0) ± |1)), to the 
final states |I/,0) and 0) respectively, at the measurement point. Since we know 
an adiabatic evolution would preserve the amplitude of probability of each state in the 
superposition, we have to design a non- adiabatic pulse in order to implement this +/- 
measurement. Indeed, by passing at an appropriate rate through the portal region, we 
achieve the desired transformation. The second plot of Fig. [6] presents the proposed 
gate. It consists of two distinct parts: the first occurs up to t < 7ns, when an adiabatic 
$c pulse is applied to the qubit. This part is done slowly to minimize the leakage when 
the system passes through the avoided crossing gap; it is also tuned so that the correct 
relative phase is accumulated between the |0) and |1) states, with respect to a reference 
phase. The second part of the pulse, t > 7ns, causes the qubit to undergo a non- 
adiabatic evolution through the portal, producing the transitions needed to implement 
the gate. The bias fiux is maintained constant during the whole gate, and it is held as 
close as is practical to the S line - in our calculation we take e = 30/i$o- The matrix 
representation Eq. ^ reveals the ideal transformation, which ends with an irrelevant 
relative phase between the ground and first excited states. For the path proposed, 
the probability of leakage is found to be 0.03%. Because the qubit has to stay, for an 
appreciable amount of time in a region of very fast relative phase accumulation, the 
main type of noise is a phase noise. The expected total net fidelity for this gate is: 
/ dS^cdSefi{S^c)KSe)F{S^c,Se) - 99.8%. 
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Figure 6. Measurement gates. (a) The pulse path used to perform 0/1 
measurement. The gate is designed to bring the system from the parking regime to 
the measurement point through an adiabatic process. The apphed flux e is maintained 
constant during the gate, and distant from the S line: e = 500/i$o- (b) The +/— 
measurement gate. The pulse evolves the system from the parking regime to the 
measurement point, performing the unitary transformation |+) = -^(|0) + |1)) \R) 
and | — ) = -^(|0) — |1)) \L). The gate consists of two different regimes: the 
first, t < 7ns, is designed to be an adiabatic process, thus minimizing leakage when 
passing through the avoided crossing. The other, t > 7ns, evolves the system in a 
non-adiabatic process through the portal, leading to the desired transitions between 
|0) and |1) states. The applied flux e is maintained constant during the gate, and close 
to the S line: e = 30/i<l>o. Insets: The gate fidelity as a function of unwanted shifts in 
<l>c and e from the optimal point of operation. The main source for loss of fidelity for 
the 0/1 measurement is related with 0/1 transitions, and for the +/— measurement it 
is related to the phase noise. Note the difference in the scales. 



4-2. Phase Gate 

The phase gate, Fig. [7| is the simplest operation of our set. In order to accumulate 
the desired phase, i. e. 9z{^t) = ujQi{t)dt^ it is sufficient to adiabatically bring 
the system from the parking point to a position where the ujqi frequency deviates few 
hundreds MHz from the reference frame; then one just has to wait the appropriate 
amount of time, corresponding to the desired phase accumulation (see inset of Fig. 
[7|. Since this position can be reached without passing through the avoided crossing, 
the leakage probability is extremely low; we calculate a leakage probability of 10"''%. 



Efficient one- and two-qubit pulsed gates for an oscillator stabilized Josephson qubit 18 




Figure 7. The phase gate, (a) The gate is designed to perform an adiabatic process, 
evolving the system from the parking regime to a position where the phase rate cjqi 
becomes appreciable compared with the reference. As one can see, that position can 
be reached above the avoided crossing gap, which greatly minimizes leakage for this 
operation. After the appropriate amount of time At at the bottom of the pulse (inset 
presents co^Oz as a function of At), the system is brought back to the initial position. 
The flux e is kept constant during the operation, (b) The fidelity of the gate as a 
function of unwanted shifts in $c and e from the optimal point of operation. The main 
source for the loss of fidelity is related with phase noise. 



Since we remain essentially within a very stable parking regime, the gate, for any value 
of is very insensitivity to fluctuations of the both applied fluxes. Thus, the gate 
can be performed with a very high total net fldelity: J dS^cdSe iJi{5^c)l^{5^)F{5^c) 5^) ^ 
99.999%. Since the operation is designed to perform an adiabatic evolution of the 
system, the 0/1 transitions are almost completely suppressed. Consequently, the main 
source of (small) fldelity loss can be characterized as a phase noise. Indeed, if we take 



a look at the operator-sum representation of the system superoperator, Eq. (12), we 
obtain the following decomposition: £"0 ^ 0.99999(7^, Ei ^ 0.0031, £"2 = and E^ = 0, 
which clearly indicates that we have just one form of noise. 

4.3. Hadamard Gate 

Completing our universal set of one-qubit gates is the Hadamard gate. Fig. [8| 
Its implementation is more complex than the previous gates, since it requires the 
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Figure 8. The Hadamard gate, (a) The gate consists of a composition of pulses in 
and e fluxes. The pulse is designed to bring the system from the memory regime to 
the portal in an adiabatic process. When the portal is reached, the system is kicked 
from one side of the S Line to the other through a non-adiabatic e pulse, in order 
to obtain the wanted transitions. After applying the e pulse, the system is brought 
adiabatically back to the memory position, (b) The fidelity of the gate as a function 
of unwanted shifts in and the synchronization of the pulses, 8t^ from the optimal 

point of operation. The main source for the loss of fidelity is related with noise phase. 



synchronization of the two parameters $c and e. As one can see, the gate can be 
basically divided in three parts: the first, t < 10ns, evolves the qubit in an adiabatic 
process from the memory state to the portal. During this process, only control fiux is 
changed; once reached the start of the portal, the second part of the gate is applied. It 
consists of a non-adiabatic pulse performed by the bias fiux e. This pulse is responsible 
to very quickly move the qubit from on side of the S hne to another. This process, 
10ns ^ t < 14ns, is designed to create the desired superpositions between states |0) 
and |1). In order to require reasonable rise-times for this pulse, we have to be as far as 
possible from the S hne. Nevertheless, this excursion in the bias fiux e is limited due to 
the very small gap region seen in the fiux space, see Fig. [5| Finally, the last step, 
t ^ 14ns, adiabatically brings the qubit back to the its initial position - again, only $c 
is changed in this regime. 

As highlighted in Fig. [5| the path taken in the fiux space due to this fiuxes $c(^) 
and e{t) composition has the form of a "U"-shape. This is a very useful feature, since 
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working on both sides of the S hne, one can expect to correct in the second half of the 
pulse some of the errors accumulated during the first. Indeed, because of the symmetry 
around the S line, in the first order of approximation, the errors due to the shifts in the 
bias fiux, 5e^ should be canceled out, since at this order 6Joi(±e + 8e) = cJoi(|e|) ± a5e. 
However, for this to be true, a precise synchronization of the $c(i) and e{t) pulses is 
required, in order to obtain a symmetric path. For the paths in Fig. [5} the number 
of dots indicates the rate as the system evolves: The fewer dots, the faster the time 
flux rate, one can clearly see the other strategies used to optimize the gate: at the 
parking regime, we can perform a very fast evolution, since the cjqi and are very 
big 3.1GHz); once we have reached the avoided crossing position, in order to avoid 
leakage, we slow down the evolution. However, this regime also coincides with a fast 
phase accumulation rate, so we try to not spend too much time at this region. Thus, 
we find a trade-off that has to be worked out during the optimization of the gate. For 
the path shown, we expect to observe a probability of leakage of 0.12%. 

The second plot of Fig. [8] shows the fidelity as a function of 5^c^ and the fiux 
desynchronization 5t. It turns out the gate proposed is very insensitive to fiuctuations 
in 5$c- This occurs because we have designed the gate to explore a "sweet" spot in 
control fiux. As one can see in Fig. [2| the region close to $c = 1.447$o presents a 
first-order insensitive point to fiuctuations of $c for several values of e. Thus, since 
this point is at the portal, we have chosen it to "sit" the qubit while we perform 
the bias pulse. The average fidelity considering each noise channel separately is 
presented in Table [T} Because the 5e and St fiuctuations essentially have the same 
effect, changing the relative phase of the two computational basis states, their loss of 
fidelity are found to be very similar. Thus, we see that the phase noise again plays 
the major role for the loss of fidelity. The complete noise characterization is obtained 
through the operator-sum representation of the system superoperator: £"0 ^ 0.9953H, 
El ^ 0.063^2^, E2 ^ 0.0048(d + {a, - (t^)/4) and E^ ^ 0.0016(d/2 + a, - a^). 
The net fidelity considering the effects of all noise channels together is given by: 
/ dS^cdSedStfi{S^c)KSe)fi{St)F{S^c,Se,St) - 99.46%. 

5. Two-qubit System 

The two-qubit system we have used to simulate the two-qubit gate is sketched in Fig. [T} 
This layout preserves the same structures present for the one-qubit system, i. e. the 
transmission line, the readout SQUIDs and the fiux lines; in addition, it has a qubit- 
qubit interaction that is assumed to arise due to the mutual inductance between the two 
big loops [H] (other qubit-qubit coupling implementations using tunable interactions 
are demonstrated in [22[ [23]). The qubit-qubit mutual inductance is considered to be 
a small parameter, such that a first-order perturbation theory is expected to give a 
very fair description of the system dynamics. Thus, following the procedure adopted to 
derive Eq. ([5]), we obtained the two-qubit system Hamiltonian given by (see Appendix) 

n =Ha + Hb + Hj, (13) 
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Figure 9. (a) The effective qubit-qubit interaction J as a function of the appUed 
fluxes ^c}' Even though the physical inductive interaction between the qubits is 
assumed constant, the effective qubit-qubit interaction is tunable as a function of the 
applied control flux in each qubit. This happens because of the change of the character 
of the system state when one changes the applied control flux in the qubit. (b) Solid 
curve presents J when one of the qubits is maintained deeply in the single- well regime. 
Dashed curve shows the case when the qubits are moved in resonance through the 
control flux space. The bias flux B' (dot-dashed curve) due to the circulating current 
of one of the qubits is presented as a function of its control flux. 



Hi = \B\<^^,mt)^^ 1^ - )1^ ® af + J($^, $f ® af . (14) 

Ha,b represent the single qubit Hamiltonians, Eq. ([5]), of the qubits A and B, 
respectively. The qubits are assumed to have identical bare qubits, but with different 
fundamental-mode transmission lines, uj^ ^ uj^. This choice was made in order to avoid 
possible double excitation of the transmission lines due to the transfer of one quantum of 
energy from one transmission line to another. In addition, because we would like to use 
the same master clock to track the dynamics of both qubits, we had to impose a rational 
relation between uo^ and uo^. In our calculations, we have assumed that uo^/uo^ = 3/4. 

The system interaction Hamiltonian Hi is given by Eq. (14). Its terms may be 
understood as arising due to "classical" and "quantum" magnetic field components, in 
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the following sense: For the first two terms, we clearly can identify them having a form 
of a bias field B' applied to each qubit (compare with the bias term of Eq. Q). As 
expected, the bias field B' applied to one qubit, let us say qubit depends of the control 
parameter of the other, qubit j. This happens because the external circulating current 
of qubit j creates a magnetic field, B'{^{)^ seen by qubit i. Nevertheless, for a given 
the qubit i feels the same bias field B'{^{) whatever the quantum state of qubit j 
is. Thus, B' should be seen as a simply result of the Faraday's law applied to classical 
circuit: The field induces a persistent current in the circuit j, which, in turn, creates 
a magnetic field B'{^{) seen by the other qubit. Not surprisingly, the first two terms of 
Hi cannot generate qubit-qubit entanglement. The dot-dashed curve in Fig. ^ shows 
B' a function of the control fiux. Because of the presence of the Josephson junctions 
(non-linear inductances), B' has non-linear behavior for some values of $c (double- well 
regime) . 

As a consequence of the existence of the magnetic field B' ^ for the two-qubit system, 
the S line of qubit i is shifted to the lines = /c$o ± ^X^c) (where k is an integer, and 
the choice of the sign depends on which loop the fiux B^ is threaded through). Thus, 
we refer to those "new" qubit symmetric lines as the S' lines. This is a very important 
eflFect and must be taken into account when performing one-qubit operations for the 
two-qubit system. 

The last term of Eq. (14) describes a magnetic field seen by one qubit, determined 
by the quantum state of the other. The result is a qubit-qubit coupling of the form 
® (rf , which is responsible for generating qubit-qubit entanglement in the system. 
As one can see in Fig. [9], even though the physical qubit-qubit inductive coupling is 
assumed fixed, they effective interaction J is tunable as a function of both control fiux 
parameters. In fact, as shown in the second plot of Fig. |9} the coupling J never reaches 
values higher than 450MHz when one qubit is maintained deeply in the single well regime 
(solid curve), rather than 5GHz observed when both qubits are at the measurement 
point (dashed curve). This is a direct manifestation of the change of the system state 
character as a function of control parameter. 

In addition, because the transmission line of qubit A sees the qubit B only via 
interaction with its bare qubit, and the interaction between the bare qubits is in fact 
small, it turns out the transmission line states are very stable under the changes of 
control parameters of the other qubit. Indeed, since the qubit-qubit interaction Eq. 
(14) is capable neither of moving the qubits from the single- to the double- well potential 
regime (for that terms of the form are needed) nor of changing the quantum number 
of the transmission lines, once one qubit is parked, the two-qubit system eigenstates 
stay "frozen" at the subspace in which the parked qubit eigenstates are close to the 
states 0) and I^S, 1). Thus, parking one qubit imposes a selection rule on the system 
that leads to a further reduction, below the already small value of J, of the effective 
qubit-qubit coupling. In fact, our calculations of the interaction matrix elements 

j^l J($^, $f )(7^ ® af|/^,m^)| (where the states iJJ.m are the 0-1 one-qubit 
states) show they never reach values higher than 3MHz when one of the qubits is parked. 
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Table 2. The expected one-qubit gate fidelities, as a function of unwanted shifts 
in 5^ci and 5t^ when performed to the two-qubit system using the same schemes 
described in Section [4j We observe a small difference between the one- and two qubit 
system fidelities only for the +/- measurement gate. 



Two-qubit system 


Gate 






8e 




0/1 Measurement 


(Fidelity) (%) 


99.99 


99.99 




+/- Measurement 


(Fidelity) (%) 


99.76 


99.79 




Phase Gate 


(Fidelity) (%) 


99.999 


99.999 




Hadamard 


(Fidelity) (%) 


99.87 


99.47 


99.63 



As a result, parking leads to a very effective way to decouple the qubit s, thus 
allowing one to perform one-qubit gates for the two-qubit system discussed. By simply 
taking account of the fact that the qubit S lines in the flux space are moved to the 
S' lines, we can use exactly the same shaped pulse schemes presented previously for 
the one-qubit system, in order to perform the universal set of one-qubit gates. Indeed, 
our simulations showed that one can expect the same probability of leakage previously 
reported, and virtually the same expected fldelities, see Table [2} 



5.1. Controlled- Z Gate 



The two-qubit gate proposed is a gate in the equivalence class ^24j of the controUed-Z 
gate, i. e. it only diflFers from the controUed-Z gate by local operations to the qubits 
A and B. The designed gate, flrst plot of Fig. [TOj involves an adiabatic evolution of 
the system states, in order to accumulate the correct relative phases. Both $c and e of 
qubits A and B are used to perform the operation. In order to obtain the strongest qubit- 
qubit interaction, and thus the shortest gate possible, the control fluxes of both qubits 
are changed simultaneously, = $f (i)- Also the bias fluxes change identically 

and, as for the Hadamard scheme, we work on both sides of the S' hue. The desired 
unitary transformation has the matrix representation in the two-qubit eigenstate basis 
{|0),|1),|2),|3)} given by 
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J 



(15) 



The relative phase 
dzz = do - 



61 — 62 + 6s 



is the pertinent parameter of this gate. Since the local invariants 
are determined by 



Gi = 1(1 + cos 6 zz)^ G2 



2{1 + ^co8 6zz). 



(16) 

of the gate Uzz 
(17) 



Efficient one- and two-qubit pulsed gates for an oscillator stabilized Josephson qubit 24 



and those of the controUed-Z gate are given by Gi = and G2 = 1, we end with the 
specific condition 6zz = 7t{1 + 2k) for the relative phase. 

Since the evolution is done adiabatically, we can write the relative phase as 
Ozzi^t) = f^^dt{{Eo{t) + E3{t)) - {Ei{t) + E2{t))}, where Ei is the instantaneous 
eigenenergy of the state i. Thus, it is evident that 6zz only has non-negligible values 



when the term Jaf ® of Eq. (14) becomes appreciable. However, as shown in 
Fig. [T0| to reach that region we have to pass through a regime of a very small gap 
300MHz) between the |3) and |4) states (minimum gap between the computational 
basis {|0), |2), |3)} and the rest of spectrum of the system). Because of that, our 
shaped dc gate has a total duration of 60ns in order to avoid leakage during the evolution. 
In addition, in order to explore "sweet" spots of operations, the designed $c(^) pulse is 
not completely flat when passing through the minimum gap region (see Fig. [To|(b) and 
(c)). 

The gate is performed in the following way: starting from the memory state, both 
control fluxes are changed to the region around $c ^ 1.453$o- Then, by changing the 
bias fluxes e^'^, the qubits are slowly moved through the region of minimum gap CJ34, 
passing from on side of the S' line to the other. At this stage, the parameter 9zz becomes 
appreciable (hundreds of MHz). The gate is designed to spend the right amount of time, 
such that the flnal evolution will satisfy the necessary condition for 9zz- Once the other 
side of the S' line is reached, the qubits are brought back to the memory state. As 
observed for the Hadamard gate, the two-qubit gate proposed also has a "U" -shape in 
the flux space. 

The probability of leakage due to this process is ^ 0.04%. Since the gaps between 
the computational states are bigger than 700MHz, unwanted transitions are strongly 
suppressed. Therefore, phase noise is the main source of loss of fldelity. As for one-bit 
gates, this phase noise arises from physical 1/f noise and pulse timing variations; we 
assume their magnitudes to be the same as for the one-qubit case, with no correlations 
between the two qubits. We then calculate (assuming, again, Gaussian noise statistics) 
the operator-sum representation of the system superoperator: 



^0 ^ V0.9905(aoi^ ® 1^ + a*a^ ® af 

+CoF®af -c*a,^®i^), (18) 



El ^ V0.0022(aii^ ® 1^ + iala^ ® a 



+CiV (8) (jf - iciaf (8) i^), (19) 



E2 ^ V0.00014(c2i^ ® + ic2a^ ® i^), 



(20) 




Figure 10. (a) ^c{t) and e{t) for the controlled-Z gate. The operation is performed 
adiabaticahy assuming ^^{t) = (t) and e^{t) = e^{t). (b) The energy sphtting 
and (c) the relative Ozz for the flux space {^^ = $f , = e^}. The energy splitting 
CJ34 gives a limit on the rate at which the system can evolve without appreciable 
leakage. The 6zz plot gives the phase accumulation rate of the principal parameter 
of the gate, Eq. (16). The path (a) taken in the flux space is illustrated with dots in 
(b) and (c). The number of dots indicates the time spent when passing through that 
region. 



where ao ^ -(1 + i)/2, cq ^ -0.1 - 0.03i, ai ^ OAi - 0.004, Ci ^ 0.6 - 0.03i and 
C2 ^ {I + 0/2 f§J The other operation elements Ek^ = 3, • • • , 16 are negligible. Thus, 



as one can see from Eqs. (18 20) these leading operation elements all lie at the subspace 

revealing the phase character of the noise. 
The total net fidelity considering the effects of all noise channels is given by 
/ dS$cdSed6tfi{6$c)KSe)fi{Si)F{S$c^Se,S^ - 99.67% 



The ideal gate for the shaped pulse Fig. 10 is such that ^0 = + ^)^^ 1^ + ^(1 - i)^^ 



, and Ej, = for k 0. 
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6. Conclusion 

Obviously, it is hoped that the results of this paper will be of direct relevance for 
experiments on qubits being performed at IBM, as well as being of general relevance 
to the problem of precision quantum control in any flux qubit system. The noise 
mechanisms analyzed here — magnetic 1/f noise, Johnson noise of circuit resistances, 
timing errors in pulse channels — are the mixture of fundamental and practical sources 
of error that are currently understood in the laboratory. Thus, the quantitative fact 
that the values of the gate infidelity (one minus the fidelity) are at the 1% level and 
below, is the major result of this paper. 

One might ask, given that this paper's main claim to importance is quantitative, 
whether the approximations made in the quantum model are actually adequate to give 
sub-1% accuracy. Perhaps the fundamental circuit Hamiltonian may be accepted as 
being very accurate; but can it be that the sequence of subsequent approximations, 
namely 1) the representation of the transmission line by a single oscillator mode, 2) the 
dimension-reduction resulting from the Born-Oppenheimer approximation, and 3) the 
final spectral truncation, involving an interpolated oscillator-displacement parameter, 
resulting in our four-level model, have the desired sub-1% accuracy? We can see, for 
example, in Fig. 4, that the spectral truncation results in changes in the absolute value 
of the energy eigenvalues, in some places by as much as 7%. 

Despite this, we believe that our infidelity numbers are nevertheless quite sound, so 
that if we calculate an infidelity here of 1%, it may be, in a more accurate calculation, 
actually equal to 0.9% or 1.1%, but not much different from that. We have a few reasons 
for saying this. First, for our gate designs, it is generally not necessary that an energy 
gap be exactly a certain value at a particular point on the flux axis; it is more important 
that, somewhere in the near vicinity of a particular flux value, the energy gap has a 
certain size. This will be true even for a model Hamiltonian of moderate accuracy. 

Traced to a deeper level, the successful functioning of our gates depends primarily on 
the following general features of the model: 1) Since many of the evolutions we consider 
are adiabatic, it is important merely that some integrated properties of the energy 
eigenvalues over some paths in parameter space be correct. This is fairly easily achieved 
by a model of moderate accuracy. 2) In a few cases basis-changing, non-adiabatic 
evolutions are important. These are achieved in just three ways: by passing through 
the portal, by crossing the S line, and (unintentionally) by moving the system in and out 
of parking. Our model obviously contains all of these features, with trends in the size 
of gaps that certainly track those of an exact calculation very closely. Furthermore, the 
general structure of the low-lying eigenvectors are captured in our truncation, meaning 
that trends in the matrix elements that determine the magnitude of Landau-Zener 
tunneling effects, are also well represented. 

So, we conclude that our qubit, as currently understood, should be capable of gate 
operation at the 1% noise level. Unfortunately, in the lab, in any given day (or month) 
there are "bugs" in the experiment that causes the qubit, often for unexplained reasons. 
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to have fidelities much worse than what we estimate here. But we nevertheless hope 
that these calculations will have significant practical value. Our gate set is universal, in 
one of two diflFerent ways, in fact: the 0/1 measurement /preparation gates, Hadamard, 
controUed-Z, and one-qubit phase gates form one well known universal set [20]; it is 
less well known that, with more overhead, Hadamard can be replaced by the +/- 
measurement /preparation gates [25j. 

So, can a "debugged" IBM qubit be used soon for universal quantum computation? 
The answer is, in our opinion, ultimately yes. The answer would certainly be no if the 
noise threshold for fault-tolerant quantum computation were in the neighborhood of the 
oft-quoted value of 10~^ |2S, EZl EBj. It is not inconceivable for the experiment to get 
to these values someday, since we find that the infidelities decrease much faster than 
linearly with the assumed noise levels. (To get to 10~^ we would need to get to the 
very daunting levels of 100n$o at IHz for the 1/f noise amplitudes and lOO/s for timing 
accuracies; there is optimism that both of these numbers are ultimately attainable.) 
Fortunately, while 10~^ was the threshold as it was understood ten years ago, much 
recent work shows that with good designs, much higher thresholds are possible [29l [30] . 
According to [30j, 1% is in fact on the high end of the noise levels for which fault 
tolerance may be possible. 

But the final answer to the question of whether fault tolerance is possible will 
require resolving a number of further uncertainties. One effect not included here is 
the fact that noise will be correlated between qubits that are physically nearby, due 
to ordinary electrical crosstalk, say. This is certainly deleterious for noise thresholds. 
However, a great virtue of our parking scheme is that, while parked, qubits are quite 
immune to electrical noise from nearby sources. Another point that might count in our 
favor is that, as seen in our operator-sum expansions, the noise has a definite structure 
(i.e., it is mostly phase noise in most cases), while most analyses of noise thresholds 
have assumed worst-case, structureless noise. On the other hand, these special noise 
sources, to the extent that they arise from the physical 1/f noise, are substantially 
correlated in time. While fault tolerance is known to survive in the presence of such 
"non-Markovian" sources [3T], it is not clear what such an eflFect does to the numerical 
values of the threshold. Finally, there are purely "architectural" considerations, e.g., 
the theoretical analyses assumed that a qubit can be moved to the proximity of any 
other one without cost. This will surely not be true in any real Josephson architecture. 

To summarize, we have demonstrated that using only low-bandwidth electrical 
pulses, a universal set of high fidelity quantum gates can be achieved for an oscillator- 
stabilized Josephson qubit, with a required "clock time" around 60ns, and achievable 
fidelities above 99% per gate operation. Essential to these results are the existence 
of "parking" made possible by coupling of the fiux qubit to a transmission line, the 
possibility of mostly adiabatic control, and the availability of two robust basis-changing 
effects obtained by crossing the "portal" and the "S line" . Time will tell if these physical 
elements make it possible, in the face of the many sources of noise that are present in 
the solid state environment, to do large-scale quantum computation. 
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Appendix A. Hamiltonian derivation 

In this appendix we supply derivations of the Hamiltonians presented in the main text, 
Eq. g and Eq. Q. The main assumption in the derivation is the use of first-order 



perturbation theory to treat terms of the system potential. As presented in Fig. [5] and 
previously discussed, performing the procedure described here, we could arrive at a very 
good description for the lowest states of the IBM qubit. 

Appendix A.l. Brief review of BKD 

Burkard, Koch and DiVincenzo [16j have introduced a universal method for analyzing 
any electrical circuit containing Josephson junctions, provided only that its elements 
can be represented by lumped elements. The methodology can be summarized in a few 
steps: A network graph is written for the circuit. A network graph is simply a drawing 
of the circuit where each two-terminal element (inductor, capacitor, Josephson junction, 
etc.) is represented as an oriented labeled branch connecting two nodes. Then, a tree 
of the network graph (a subgraph that does not contain any loops) is chosen using the 
following criteria: the tree has to contain all of the capacitors in the system, no resistors 
or external impedances, no current sources, no Josephson junctions and as few linear 
inductors as possible. The branches do not belong to the tree are called chords. 

Associated with the tree chosen, there are the so-called sub-loop matrices F^y, 
where the label X represents the tree elements {X = C for capacitors and X = K ioi 
the tree inductors) and Y the chord branches (Y = J for Josephson junctions, Y = L 
for linear inductors, Y = R ioi shunt resistors, Y = Z ioi external impedances and 
Y = B for bias current sources). The sub-loop matrices have entries -1, 0, or 1, and give 
information about the interconnections in the circuit, determining which tree branches 
X are present in which loop defined by the chords Y . The entries in F^y matrix are 
found as follows: If the corresponding tree element Xi does not belongs to the loop 
defined by the corresponding chord Y^-, its entry F^y is zero. If it belongs to the loop 
but has opposite orientation to 1^, we have Fj^y = +1. Finally, if it belongs and has 
the same orientation, we have F^y = — 1. These steps give an algorithm to encode the 
topology of the system in a matrix representation. In addition, the formalism assumes 
that all capacitors should be considered to be in parallel with a Josephson junction, 
even if it is one with zero critical current. 

The physics of the circuit is introduced by imposing the KirchhoflF laws at each node 
of the network, and defining the electrical characteristics of each branch type as 

Ij =lcsin(^, (A.l) 

Qc = CVc, (A.2) 

= RI^, (A.3) 

Yz{uj) =Z{uj)lz{uj), (A.4) 
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T ^ 




(A.5) 



Here the diagonal matrix Ic contains the critical currents of the junctions, and sincp 
is the vector (sin 9^1, sin 9^2, • • • , sin (/pa^j). The phase (fi represents the superconducting 



phase difference across the junction i. The linear capacitors are described by Eq. (A.2), 



with their capacitance values given by the entries of the diagonal matrix C. The junction 



resistors are assumed to follow Ohm's law, Eq. (A.3), where R is the real and diagonal 



shunt resistance matrix. The impedances are described by Eq. (A.4), that relates the 



Fourier transforms of the current and voltage. Z(cj) is the diagonal impedance matrix. 
Since linear inductors can be tree branches as well as chords, we have to distinguish 
them to apply the network graph theory. Therefore, the inductance matrix must be 



organized in the block form shown in Eq. (A.5), where L is the inductance matrix of 



the tree inductors, is that for the chord inductors, and L^k represents the mutual 
inductances between the tree and chord inductors. 
BKD arrived at the system Hamiltonian 

l-CfcC-'Qc+(p-Xu{^,t), (A.6) 



nsit) 

U{cp,t) 



1 271 

J2 L-ji cos + -v?^MoV? + ^^''^[(N * + (S * lum], (A.7) 



where the external applie d flu xes and current are respectively represented by and 
Ib- The first term of Eq. (A.7) is due to the Josephson energies, where Lj^ = I^-i 
and /c;i is the critical current of junction i. The second is associated with the inductive 
energies of each branch of the circuit. The topology of the circuit is encoded in the 
matrices matrices Mq, N and S, given by Eqs. (62), (63) and (66) in [16]. Thus, BKD 
maps the circuit dynamics to that of a massive particle in a potential, whose masses and 
degrees of freedom are associated with the system capacitances. The quantization of the 
system is introduced by imposing the canonical commutation relation for the variables 
of charge, Qc, and phase cp: Qc-j] = ihSij. 



Appendix A.2. IBM qubit network graph 

Fig. |A1| shows the network graph and the chosen tree for the bare qubit. The sub-loop 
matrices can be read off by inspection: 
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Figure Al. The network graph for the bare qubit. Each Josephson junction Ji 
(thick hue) is modeled by a resistively shunted Josephson junction (RSJ) containing 
an ideal junction with critical current /c,i, a parallel junction capacitance Ci and a 
shunt resistance Ri. The qubit is coupled to the flux lines via mutual inductances M12 
and M34. Because of the "gradiometer" structure of the qubit, the "pick up" flux line 
(responsible to apply the flux $p) was omitted without loss of generality. Inset (b): 
tree chosen for the graph. 



The L matrices are denoted: 
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(A.12) 



Using the modeling program FASTHENRY, we obtained the foUowing inductance 
values for our qubit: Li = 32.18, L3 = L5 = 686.93, L2 = 32.18, L4 = 605.03, = 
106.27, Mi3 = Mi5 = -2.23, M35 = -12.25, Mu = 0.8, M34 = 0.5, Mi^ = 
O(exact), M3C = -M5C = 28.15 (all in units of pH). 

In addition, the capacitances and critical currents of the junctions are taken to be 
C = lOfF and = 1.3M- 
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Appendix A. 3. Bare qubit Hamiltonian 

Our aim in this section is to present the steps used to derive the approximate 
Hamiltonian for the bare qubit. In order to reach the appropriate conditions to perform 
a first-order perturbation theory to treat the couphng between the qubit and the bias 
fiux hne, we firstly have to identify the "slow" and "fast" coordinates of the system. 
The slow coordinate is that that joins the two minima when the system presents a 
double-well structure, or that which has the slowest curvature in the single- well regime. 
The fast coordinates are those in which the potential rises very steeply, such that the 
system dynamics is frozen into the ground state along these directions. 

The procedure to identify the slow and fast degrees of freedom starts from the exact 
system potential Eq. ^ and a unitary transformation R that diagonalizes Mq (Mq 
is a real and symmetric matrix [16j). Using R we can decouple the quadratic form of 
the potential by transforming the system coordinate to the new coordinates n = Rep. 
Thus, from Eq. ^ we obtain (ignoring the term due to external sources of current) 



[/(n, t) = - ^ LjI cos rikU]^ 



+ VAn+(^)n^$., (A.13) 

where we have defined the diagonal matrix A = RMqR^, which entries are the exact 

eigenvalues Ai = 0, A2 = 3/(L3 + 21^ - AM^^ - M35) and A3 = (Li + 2{L^ + M35 - 

2Mi5))/(Li(L3 + M35) - 2Ml^), and the "non-orthogonal" flux vector = RN*^; = 
f A A2 , ^(i:^3-Mi5+M35) ^ -1 -^X^ 

i^' 75^' -Li(L3+M35)-2M?5'^c|, Wltil 

e =$-$p, (A.14) 

Because of 2(L3^M3^^Mi5) ^ 0.024 and the fact we expect to work with magnetic fiuxes 
of the order of only few $oS, we have that, to good approximation, the system potential 
is only a function of the diflFerence of the magnetic fiux in the two large loops, defined 
as e. This property is a direct manifestation of fact the designed bare qubit has a 
"gradiometer" structure. Even though this correction is expected be very small, we 
shall consider it during the whole procedure we are going to describe here. 

The new coordinate system allows us to easily identify, by inspection of each term 



of Eq. (A.13), the following symmetry in the potential 

[/ (n, l>e, $, $p) = [/(ni, 712,713, $c,e) 

= U (-ni, -n2, 713, -e) . (A.16) 

This shows that the system potential presents a definite symmetry in the plane 
determined by the directions {711,712}, when the system goes from +6 to — e (if we 
simultaneously adjust the control fiux $c to maintain the new fiux coordinate Eq. 
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(A. 15), constant). In addition, Eq. (A. 16) also permits us to determine the conditions 
for which the system presents the same physical potential 

U (ni, n2, 713, *c, e) = U {n\,n'^, n'^, e') + 7, (A.17) 

where 7 is a irrelevant constant. In fact, we find that whenever the bias fiux difference 
Ae = — e is an integer multiple of the fiux quantum, let us say Ae = {ki — /c2 — 2/c3)$o 
(each ki is any integer), the Eq. ( |A.17D is satisfied if we have the fiux shift A$c = 
— $c = ^^2(Lf-f Afe-Mit)^^ ^o(fei + ^2) accompanied by the Josephson phase changes 
A(/Pi = 27r/ci, A(/P2 = 27r/c2, and A993 = 27r/c3. 

It is now worth introducing a shifted coordinate system defined by n = s + d, where 
d = {0, 0, ^ (^g) l^StSf^*^^' potential written in this coordinate system 
is given by 

1 - 2 , A 2\ , f 27r\ A2 



U (s, e) = -{X^si + Xssi) + j —6^2 

Observe that the shifted coordinate system s is such that we still have the symmetry 
observed in Eq. (A. 16): [/(^i, 52, 53, $c, e) = ^^i, — 52, 53, — ^)- However, since 
A2/A3 ^ 0.059 and Lj.l/X^ ^ 0.061 {i. e. the potential is much steeper in the direction 
53 than in the others), we have the global minimum of the potential occurring at the 
position 53 0. As a result, the bare qubit low-level dynamics can be considered frozen 
at the plane 53 0. 

Another important result is derived from the fact the low-level system dynamics is 
only governed by the degrees of freedom Si and 52 : the system potential has a symmetry 
line (S line) at e = in the fiux space. Indeed, as one can see, when e = 0, [/(s, ^) 
has a perfect symmetry around the origin in the {51,52} plane. Because of that, the 
low-level system wave functions have definite parity symmetry at e = 0. In addition, 
because the same physical potential is found when e is changed by an integer multiple 
of $0, a S line occurs whenever e is an integer multiple of $o- 

Unfortunately, because A2 ^ L'J^^ it is not simple, as it was for the direction 53, to 
determine the soft and fast directions in the plane defined by 53 0. However, from 
the symmetry of the S line one can see that, if the potential has a minimum in 5^^^ 7^ 
and 5^^^ 7^ 0, the potential presents a symmetric double- well structure (with maximum 
at 5i = 52 = 0) in the direction connecting the minimum points. This direction is 
expected to be the slow coordinate of the system, while its perpendicular direction should 
determine the fast degree of freedom. Consequently, by finding the minima positions 
in the plane {^i, 52}, one can determine those directions. Following this procedure, we 
perform one more rotation, qi = Si cos 9 — S2 sin0, q2 = Si sinO + 52 cos 9 and qs = 53, 
such that the direction qi connects the minima. When the potential does not have a 
double- well structure, we define the last rotation so that the potential curvature in the 
direction qi is the smallest. 
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Thus, after appropriate transformations, we end with phase and flux coordinate 
systems in which the symmetries presented by the system are much more clearly stated. 
In addition, the slow, qi and fast, {^2,^3}, degrees of freedom of the system are also 
identifled. As already mentioned, because the potential is very steep the fast directions, 
the low-level system dynamics is frozen into the ground state along these direction. 
Following the Born-Oppenheimer approach developed in [13], those directions can be 
traced out and their eflFects incorporated as small corrections to the remaining slow- 
coordinate potential energy. 

At last, following all the steps described above, we are now in position to construct 
term by term of Hamiltonian Eq. ([5]): 

• No bias flux line (e = 0): As the system potential has a perfect symmetry around 
the origin, the bare qubit wave functions also have a parity symmetry. Thus, the 
ground and flrst excited states are expected to be symmetric and antisymmetric 
with respect to reflection around the origin. As described in the main text, the 
system potential is found in a double- well structure for small values of $c- There 
the bare qubit eigenstates can be understood as symmetrical and antisymmetrical 
equal superpositions of the classical orbital states of the left and right wells. For 
the single- well regime, $c ^ l-S^o, the bare qubit states are associated with those 
of a harmonic oscillator. Thus, in this representation, the Hamiltonian of the bare 
qubit without bias flux line can be written as —\/\{^c)^x) where A($c) (shown in 
Fig. [3]) is the ground and flrst excited state energy splitting. 

• The bias flux term: As previously discussed, the bias flux e is responsible for 
breaking the perfect potential symmetry observed at the S line. However, the 
system presents, as stated in Eq. (A. 16), another important symmetry in going 
from one side of the S line to another [i. e. passing from —e to +6): the 
system potential at e < is identical to that in e > by a mirror reflection (a 
TT rotation for the full multi-dimensional potential) at the origin. In the double- 
well regime that corresponds to the interchange of the left and right states under 
the mentioned transformation. The bias term of the Hamiltonian Eq. ^ is 
obtained applying a flrst-order perturbation theory to treat the bare qubit-bias 

flux line coupling. Indeed, because the potential term representing their coupling, 
2 



Uq^ ^ egiTry IA2 sin0, is very small, for typical values of bias flux used in this work 
(e = 0(m$o))5 compared to the other potential terms, the flrst-order approximation 
can be done in a very controllable way. Because of the symmetry of the bare 
qubit eigenstates, the interaction matrix elements (where = S^A are 

the symmetrical and antisymmetrical bare qubit eigenstates) are non-zero only for 
those connecting the ground and excite states. Thus, the bias term has the form: 
|6($c)^z, with 

6($,) = 271 sin 9 {S\q,\ A)' (A.19) 
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Appendix A. 4- Bare qubit-transmission line coupling 

Another structure present in the IBM qubit is a high-quahty superconducting 
transmission hne Fig. [TJ We model the fundamental mode of the transmission line 
as a simple LC circuit of definite frequency ujt coupled to the bare qubit. We assume 
its characteristic impedance and inductance are given by Zo = = llQVt and 

Lt = 5.6nH respectively, and the bare qubit-transmission line mutual inductance equals 
MqT = 200pH. In order to find the slow and fast coordinates, we apply the same steps 
discussed above. An interesting feature of the new system is the switch of the slowest 
coordinate from the bare qubit degree of freedom to the transmission line coordinate. 
Indeed, we observe for small values of when the bare qubit ground and first excited 
state energy splitting is much less than HljUt^ that the slow coordinate is associated with 
that of the bare qubit. However, when the energy splitting is higher than Hujt^ at high 
values of the slowest coordinate is that associated with the transmission line degree 
of freedom. Because of that, when applying the Born-Oppenheimer approximation, for 
a correct description, we stop with a two dimensional potential energy of coordinates qi 
for the bare qubit and qT for the transmission line. 

Provided that the bare qubit-transmission line coupling UqT is very small compared 
to the other potential terms: UqT/\2 ^ 10~^, we can safely apply a first-order 
perturbation theory to treat that coupling. Since the interaction term has the form 
UqT oc qiqr^ the only non-zero interaction matrix elements (where 
^ = S^A are the symmetrical and antisymmetrical bare qubit eigenstates, and Lp = 
0, 1, • • • , n are the harmonic oscillator states) are those connecting the ground and first 
excited bare qubit states and the transmission line states differing by one quantum of 
energy. Thereby, we obtain the Hamiltonian form: g{^c){^ + ^^)^z^ with 

g{^,) ^ {{S\qM)mQT\l)) (A.20) 



Appendix A. 5. Qubit-qubit coupling 

The layout of the two-qubit system is shown in Fig. [TJ This scheme contains the same 
structures as the one-qubit system: the transmission line, the readout SQUIDs and 
the fiux lines. The bare qubits are assumed identical, but, in order to avoid possible 
double excitations of the transmission line states due to the transference of quanta of 
energy from one transmission line to another, the transmission line fundamental mode 
frequencies are made different uj^ ^ uj^. As highlighted in Fig. [ij the qubit-qubit 
interaction arises due to the mutual interaction between the two big loops, which is 
assumed to be by M = 12pH. 

Following the steps introduced for the bare qubit Hamiltonian, we reach after the 
diagonalization of the Mq matrix the system potential 

U{n^, n^, t) = [/(n^, t) + [/(n^, t) + [/^^ (A.21) 

where U{n\t) corresponds to the potential Eq. (A. 13), and [/^^ = ^ijCijufn^ is 



the interaction potential. We find that c^^^/A2 ^ 10 ^, what justifies considering the 
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qubit-qubit coupling a small parameter of the potential, and consequently that the qubit 
dynamics is weakly perturbed by their interaction. Thus, we can perform the same steps 
described for the one-qubit system to find the qubit s slow and fast coordinates. 

An interesting result arises when the qubit potentials are written in the shifted 



coordinate system Eq. (A. 18). Under this transformation, the interaction potential 
^AB cj^anges to the form 

^AB 



(A.22) 



where 4 = ^ r^^SrSFt^*'^' "frozen" value of the fastest degree of 

freedom of qubit i. After performing the final rotation to align the slow degree of 
freedom with one of the coordinate directions, the pertinent terms of the interaction 
potential are given by 

U^^ = ag^^gf + /3(rff sin^V - rf^sin^^gf). (A.23) 



As one can see, the linear terms in Eq. (A.23) are due to the "frozen" values of the 



fastest degree of freedom of each qubits. Moreover, the linear term of the slow coordinate 
of qubit i arises due to the frozen coordinate of qubit j and vice- versa. The interaction 
matrix elements due to these terms are exactly that of the one-qubit bias term, leading 
to a Hamiltonian form: 



(J. 



Finally, the term QiQi of Eq. (A.23) only connects states of diflFerent parity in both 
qubits. Thus, this term leads the form: 
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